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Abstract 

In this paper we detail some results advanced in a recent letter [Phys. Rev. Lett. 102, 073008 (2009) ] showing how 
to engineer reservoirs for two-level systems at absolute zero by means of a time-dependent master equation leading to a 
nonstationary superposition equilibrium state. We also present a general recipe showing how to build nonadiabatic coherent 
evolutions of a fermionic system interacting with a bosonic mode and investigate the influence of thermal reservoirs at finite 
temperature on the fidelity of the protected superposition state. Our analytical results are supported by numerical analysis of 
the full Hamiltonian model. 



I. INTRODUCTION 

The purpose of the engineering reservoir program [T] is 
to protect a specific state against the decoherence stem- 
ming from the natural coupling between a quantum sys- 
tem and the reservoir. To engineer a reservoir, a given 
system, whose state is to be protected, is compelled to 
engage in additional interactions besides that with the 
natural reservoir. The engineering reservoir technique is 
then applied to make these additional interactions pre- 
vail, modifying the dissipative Liouvillian in a specific 
way that drives the system to the desired equilibrium 
with the engineered reservoir. Rightly connected to the 
engineering Hamiltonian program [5], the reservoir engi- 
neering has been developed for trapped ions O |3] and 
atomic two-level (TL) systems [5]- Recently, under the 
assumption of a squeezed engineered reservoir, a way to 
observe the adiabatic geometric phase acquired by a pro- 
tected state has been proposed (6n8j under the assump- 
tion of an adiabatic evolution of the system-reservoir pa- 
rameters. It is worth mentioning recent results on the 
possibility to drive an open many-body quantum system 
into a given pure state by an appropriate design of the 
system-reservoir coupling [HI HD] . 

In this paper we detail the general recipe presented in 
Ref. [llj showing how to build nonadiabatic coherent evo- 
lutions of a two-level ion trapped into a leaky cavity, in 
contrast with the applications of the engineering reservoir 
technique which require adiabatic conditions ^-'H'. Fur- 
thermore, we show how to implement the proposal of Ref. 
[T^ for controlling the slow changes of a protected TL 
system through the parameters of the engineered reser- 
voir. We also investigate the influence of thermal reser- 
voirs at finite temperature on the fidelity of the protected 
superposition state. 



This paper is organized as follows. In Section II we 
briefly present the engineering reservoir program, includ- 
ing the advances presented in Ref. [TT] . In Section III we 
show in details a general recipe to engineer reservoirs for 
a fermionic system interacting with a bosonic mode. As 
an application we derive an effective reservoir for a time- 
dependent superposition state of the internal degrees of 
freedom of an ion trapped into a leaky cavity. In Section 
IV we present numerical results supporting our approx- 
imations and an analysis of the temperature effects on 
the protected time-dependent superposition. Finally, in 
Section V we present our conclusions. 



II. THE EXTENDED ENGINEERING RESER- 
VOIR PROGRAM 

Before we present the technique to protect superpo- 
sitions of quantum states evolving nonadiabatically, we 
review basics results concerning the engineering reservoir 
program. To this end, we focus our attention on a system 
interacting with a reservoir, both system and reservoir 
being modeled by harmonic oscillators. The coupling, as 
usual, is assumed to be linear in position-position. From 
this model one can deduce a master equation, obtained 
by tracing over the reservoir variables, to study the evo- 
lution of the single harmonic oscillator. By driving the 
system with external flelds and/or allowing the system to 
interact with additional systems, each one possessing its 
own reservoir, it is possible to obtain, through some ap- 
proximations, an effective master equation which is the 
starting point to protect a given state by means of engi- 
neered reservoirs, as we detail in the following. 



A. Protection of a stationary quantum state 

The master equation describing the evolution of the 
density operator p of a given system coupled with its 
natural reservoir in the interaction picture and at zero 
temperature is 
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where (a) is the creation (annihilation) operator in 
the Fock states of the system. As it is well known, the 
only steady state resulting from Eq.Q is the vacuum |0), 
which is an eigenstate of the annihilation operator with 
zero eigenvalue: a |0) — a |0), a = 0. The main goal of 
the standard engineering reservoir [1] is to obtain, in the 
interaction picture and at zero temperature, an effective 
master equation in the form 



{20pO^ - O'^Op - pO'^O) , 



(2) 



where V is the effective decay rate of the engineered reser- 
voir which is coupled to the quantum system in a specific 
way characterized by the time-independent system oper- 
ator O. Proceeding in analogy with Eq.Q, the only pure 
steady state of this system is the eigenstate \'ip) of the op- 
erator O with null eigenvalue, ensuring that there is no 
further eigenstate \4>) of O such that [0,0^'] \4>) = 0. As 
a consequence, is the asymptotic state of the system 

a- 



engineered Hamiltonian H (t) prompts the specific oper- 
ator O {t) and its protected evolving eigenstate \^ {t)) 
(with null eigenvalue). The key feature to be noted here 
is that, through the unitary transformation R (t), we re- 
cover the time-independent form of the master equation: 

PR = ^ (20rprOI, - 0\,Orpr - pr0^rOr}i . (4) 

Therefore, the protected stationary state, {iPr) 
(OrIiPr) — 0), turns out to be a nonstationary 
state \ip(t)) — Rit) I^-r) in the original interaction pic- 
ture of y, where H (t) = iR{t)R'f{t). Getting rid of the 
adiabatic constraints, we are thus allowed to manipulate 
the evolution of the protected state \ip{t)) through appro- 
priate engineered Hamiltonian and reservoir. It is worth 
mention that the time dependence of the protected state 
is closely related to the properties of H{t) and O (t). If 
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is straightforward to see that |V'(i)) — IV'-r): i-^-, IV'(O) 
stationary apart from a global phase factor. 

Next we develop a general theory of reservoir engineer- 
ing for a system interacting with a bosonic mode, i. e., 
we show how to use the bosonic decay to build up the 
general master equation ([2]). Then we apply this theory 
to protect quantum states of a TL system, showing how 
to obtain an effective interaction between this system and 
a cavity mode which leads to the desired state protection 
mechanism. 



B. Protection of a time-dependent quantum state 

The Ref. adds an improvement on the engineering 
reservoir program by allowing to remove the adiabatic 
constraint in the decoherence-free evolution mentioned 
above. To understand this simple, yet effective improve- 
ment, here we review the theory developed in Ref. 
Consider the engineered time-dependent master equation 
in the interaction picture 

P^~i[H it) , p] + ^ [20 it) pO^ (t) - Ot [t) O it) p 
-pOt(i)0(i)] , (3) 

where the Hermitian Hamiltonian H (t) must be en- 
gineered aiming to synthesize — without adiabatic 
impositions — the desired time-dependence of the op- 
erator 0{t) = R{t)ORR^t), with Or = 0, R{t) = 
Texp^— i H {t') dt'^ , and T being the time-ordering 

operator. Such a relation between O (t) and R{t) justifies 
the above mentioned intimate connection between both 
programs of engineering Hamiltonians and reservoirs; in 
fact, the reservoir engineering technique relies on engi- 
neered Hamiltonians. We stress that we have achieved 
the time evolution of Eq. p| in a particular way that the 



III. RESERVOIR ENGINEERING FOR A 
FERMIONIC SYSTEM INTERACTING WITH 
A DISSIPATIVE BOSONIC FIELD 

The starting point to develop a general theory of reser- 
voir engineering for a fermionic system (from now on "the 
system" ) interacting with a dissipative bosonic field is the 
engineered effective Hamiltonian 

Heff ^ Kff {Oa^ +D^a) , (5) 

where (a) is the creation (annihilation) operator of 
the bosonic field, while D and are operators associ- 
ated with the system whose state we wish to protect, and 
Xeff is the effective coupling between the bosonic mode 
and the system. As an example, some of us showed how 
to build up bimodal interactions in cavity quantum elec- 
trodynamics using three- or two-level atoms |2] . We note 
that the interaction ([s]) is a bilinear form similar to the in- 
teraction of a bosonic field with a natural reservoir which 
leads to the well-know Liouvillian describing the ampli- 
tude damping mechanism. Thus, the engineered interac- 
tion ([5| must generate an effective dissipative Liouvillian 
which competes with the natural one. When the engi- 
neered decay rate is significantly larger than the natural 
one, the effective dissipative Liouvillian will govern the 
dynamics of the system. 
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Once achieved the important step of building the inter- 
action between the system of interest and the quantized 
bosonic mode, we now add the dissipative mechanism of 
the bosonic mode through the master equation 
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written in the same representation we obtained Hefj; 
the factor k stands for the natural decay rate of the 
bosonic mode and the supra index in p^*") indicates its 
dependence from both the field (a) and system (s) op- 
erators. From the above equation we straightforwardly 
derive the evolution equation for the matrix elements 
Pnm — {n\ p^^°''^ \m) in the Fock basis: 
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To engineer the reservoir, we assume that the decay con- 
stant K of the bosonic mode is significantly larger than 
the effective couphng Ag// in Eq. ([t]). This condition 
is easily achieved for cavities with low quality factor 
Q = lu/k, where uj is the bosonic mode frequency. To- 
gether with the good approximation of a reservoir at zero 
temperature, this regime enables to consider only the ma- 
trix elements pmn inside the subspace {|0) , |1)} of Fock 
states. Actually, Eq. Q can very well be described by 
the set of equations 
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p['l^-^Kf[Opl,:l^p[:lO^)~np[% (10) 

which are similar to those written in the atomic ba- 
sis in the reservoir engineering program for trapped 
ions [3]. The strong decay rate k enables the adia- 
batic elimination of the elements poi and pn from the 
equations above. The formal adiabatic elimination is 

equivalent to assume p^'^\q — 0, allowing us to write 



Pi.o 



-i2- 



- (^Opq^I — Pi.i^^ ■ We can thus ehminate 

(s) s) 

both p\ g and ( from the dynamics of the diagonal ele- 
ments of the density matrix by substituting them in the 
above equations. Since p^^l + p[^\ ~ p^*' = Psys, the 
reduced density operator for the system of interest thus 
results 



(2£)p,.ysD^ - D'fDpsys - PsysO^^O) , (11) 



where Teng = ^A^j^/k represents the effective decay rate 
corresponding to the engineered reservoir. With this re- 
sult we can easily see that the reservoir engineering tech- 
nique depends basically on the manipulation of the inter- 
action between the bosonic mode and the system whose 
state we wish to protect. 

The two kinds of master equations presented above, 
(H and may be recovered through Eq. (|11[) by 



identifying the operator D with O and Oh, respectively. 
Such correspondence is established when the operator O 
is written in the usual interaction picture, while On is 
written in an arbitrary representation where it is sta- 
tionary and Eq. ([3| remains valid. Below we show how 
to engineer a reservoir which allows to protect an arbi- 
trary evolving superposition state of the internal degrees 
of freedom of a trapped ion/atom without the adiabatic 
constraints. 



A. Engineering reservoirs for a two-level system 
in a leaky cavity 

To implement the ideas discussed in the previous Sec- 
tion, we use a two-level (TL) system (s) characterized by 
the transition frequency ujq between the ground \g) and 
excited |e) states. This system can be either a TL neu- 
tral atom in a dipole trap [13) or a TL ion in a harmonic 
trap with frequency v |14| . The transition \g) -H- |e) is 
driven by a classical field of frequency ujc, with coupling 
strength flc. The TL system is made to interact with a 
cavity mode field (a) of frequency oja under the Jaynes- 
Cummings Hamiltonian with Rabi frequency g. Assum- 
ing the case of a two-level trapped ion, the Hamiltonian 
modelling this system is given, within the rotating-wave 
approximation (RWA), by 

g cos \ k . X j aa^g 
+r],e'(^-^+^=""=*) (Teg + H.c.j , (12) 



where W (b) is the creation (annihilation) operator of 
the vibrational mode whose position operator is ^ — 
(6^ -|- b)/\/2mi'x, m being the ionic mass and x the unit 
vector along the vibrational direction. The wave vectors 
k and k c stand for the cavity mode and the classical 
amplification field (with relative phase 4>c), respectively, 
while ars = \f) {s\ (r and s labeling the states g or e). 
The vibrational mode is decoupled from the remaining 
degrees of freedom of our model by assuming the wave 
vectors k and fc c to be perpendicular to ~af ; otherwise, 
a sufficiently small Lamb-Dicke parameter keeps the mo- 
tional state almost unchanged. Under this assumption 
we arrive at the Hamiltonian 

H = uJaa^a+Y<^z + [gaaeg + e'^'^"-''-'^ a^g + H.c.} , 

(13) 

which will be the starting point for our purposes. In the 
following, we show how to engineer reservoirs suitable for 
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obtaining nonadiabatic evolutions of the internal ionic 
states. 



1. Nonadiabatic Evolution 

To engineer the appropriate interaction between the 
TL ion and the cavity mode we apply the unitary trans- 
formation given by Ui — exp^—i (uaa^a + UcO'z/'^) t], 
leading to the Hamiltonian 



are the decay rates of the cavity mode and the TL system, 
respectively. Now, except by the superoperator describ- 
ing the decay of the two-level system, we can see that 
Eq. (16) is equivalent to Eq. (rol), with Xeff — geff and 

£) — e^" . As described above, assuming a strong 

decay rate k we can adiabatically eliminate the cavity 
mode variables. This procedure is needed to make clear 
what is the asymptotically protected atomic superposi- 
tion. Then, the dynamics of the reduced density operator 
for the TL system results to be (p^*-' = 'pat) 



Hi 



iS„t 



H.C.] , 



with Ac = ujQ — LUc {6a — i^a ^ ^c) being the de- 
tuning between the atomic transition (cavity mode) 
and the laser field. Moving to another frame of ref- 
erence defined by the unitary transformation U2 — 
exp {-I [(Ac/2) a, + {a,g e"^- +age e^'^^] ' the 
foregoing calculations are significantly simplified. This 
procedure leads to the Hamiltonian 

i/2 = .ge-**»*a[A(a++ -a__) + 

A+_ c'2?* a+_ + A_+ c-'2«t ^_^] _^H.c., (14) 



where A 
and ^ 



- ^4^/4, A 



T(2±x)/4, 



ilcl +A2/4. The atomic operators 
are defined by <t±± = |±) (±|, with |±) = 
\ (y2±^ |e) ± e-»^= ./2Wx Iff)) and x = K./^- By as- 
suming a large detuning between the cavity mode and 
the laser field, such that Sa ^ g, together with the ad- 
ditional choice 5a = —2^, we obtain under the RWA the 
effective Hamiltonian 



(15) 



with geff = g {1 — x/2) /2. The effective coupling thus 
depends on the parameter Xy whose value follows from 
the laser detuning Ac which must be significantly smaller 
than the Rabi frequency ilc- 

Next, observing that x G [— 2, 2] , we focus on the values 
X ^ 2 which allow for an effective coupling of the same 
order of the atom-cavity field coupling, i.e., geff (x) ~ g- 
Once achieved the building of the interaction between 
the two-level system and the cavity mode, we now take 
into account their interaction with a thermal reservoir at 
temperature T ~ OK through the master equation 



Heff,P 
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(16) 
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where again the supra index in denotes the density 
operator for both system and field and the tilde is used 
to describe the operators in the same representation of 



Heff, i.e., O = UlulOUiU2. Here the constants k and 7 



Pat 



so that r. 



(2o-+-PaiO--+ - O- Pat - PatO- ) + Cpat, 

(17) 

, = Ag'^jjl K represents the effective damping 



of the engineered reservoir and Cpat is the last term in the 
r.h.s. of Eq. (16) after tracing over the cavity field vari- 
ables. Note that the operator cr^ and consequently the 

protected state |+) [a^ |+) = 0), exhibit no temporal 

dependences in the convenient representation where we 
have described the evolution ([l7|. However, as we shall 
see in the following, a coherent nonadiabatic evolution 
is recovered in the interaction picture, where |-|- (t)) = 
UU1U2 \+) and U = exp \i (wqa'^ a -|- a;oCz/2) t\ , even un- 
der spontaneous decay, provided that 7 ^ Teng ■ Here we 
note that the combined unitary operations UU1U2 act on 
the Hilbert space of the atom and the mode. However, 
once we are interested in the two-level system only and 
since the field and atomic operators commute with each 
other, we are omitting, for convenience, the correspond- 
ing state of the mode. 

Ideal case: In the ideal case where 7 = 0, the solution 
of Eq. ([T7| leads to the TL steady state pat = |+) (+|- 
By returning to the interaction picture we thus obtain 

Pat 

|+(t)) {+{t)\, where 



l+W)-^(v^k)+e-'*(*^V^l5)), (18) 



with $(i) — 0c — Act. This state describes a nonadiabatic 
evolution that depends on the detuning Ac between the 
atomic transition and the classical field. Note that under 
the restriction x ^ 2 imposed above, the states around 
the north pole of the Bloch sphere are not allowed steady 
states. 

Nonideal case: To appreciate the robustness of the 
present reservoir engineering technique, it is necessary to 
take into account the damping stemming from the natu- 
ral reservoir when 7 7^ 0. Here we will consider the reser- 
voir at zero temperature. The effect of finite temperature 
in the fidelity of the protected state will be analyzed in 
the next Section. For convenience, we will analyze such 
effects in the frame on which we have defined Hi . To this 
end, we use pat = UlpatU2 in Eq. (17), which can thus 
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be written as 



Pat 



),Pat\ 



^ (2cr+_PatCT-+ - (T Pat-Pat<y ) (19) 



2 i'^^gePat'^eg ~ ^^eePat ~ Pat<^ee) 



By projecting Eq. (191 on the atomic basis we find the 



following set of equations corresponding to the matrix 
elements p++ = {+\pat\+) , ptf = {+\Pat\-), and 
Pat'^ = {-\Pat 

pt^ = ^eng + 7 (x' " 2) - [T,„g + 2^ (x' + 4)] p++ 



. +- 

Pat 



-7X^4 - (ptt + Pat^) , 
7 (x' - 2) - 27X^4 - x^plt 



eng 



+ 27 - 12) - *2e 



pit 



+7^4 - X^Pat^ 



with 7 — 7/I6. Imposing the condition = 0, we 
can also determine the asymptotic solutions for and 
greater simplified considering a large cooperative 
parameter C = <?^/7k — ^eng/l ^ 1, given by 



P++ = 1 - e++ , 

P+- = 

where e++ = (7/re„g) [(2 + x) /8]^ and e+_ = j/Ageff. 
Under the condition Teng/l ^ 1, we see that both 

e++ and arc much smaller than unity. Actually, 

taking into account the noise effects introduced by the 
reservoir, the steady state is approximately described by 
Pat — 1+) (+1- The net effect of the noise introduced by 
the atomic decaying mechanism is computed through the 
fidelity 



J' = Tr 



IPatl 



1 



The approximation leading to Eq. ( 19 1 is better as 



higher the decay rate k; however, for the dynamics of 
the TL system to be driven by the engineered reservoir, 
the magnitude of k must be chosen within the restric- 
tion C ^ 1. Within the optical regime [TS], for example, 
where g ~ 7 x 10® and 7 ~ 2 x 10^ s""'^, we obtain for 
a cavity decay constant k ~ Sg, the strength k/j ~ 10 
and a fidelity around 96%. Therefore, under the excellent 
approximation e+_|_ <C 1 and reversing the unitary trans- 
formation R{t), the state \+), written in the interaction 
picture Eq. ( 18 ) allows for a nonadiabatic coherent evo- 



lution of a TL system under spontaneous decay, which 
can be manipulated through the parameter Ac. Such an 
evolution corresponds to flips in the atomic states, rep- 
resenting trajectories on different planes parallel to the 
equator on the Bloch sphere, governed by the Hamilto- 
nian H = fl^aen e-^^^-*""^-) +H.c.. 



IV. NUMERICAL RESULTS AND TEMPERA- 
TURE EFFECTS 

To validate the approximations carried out in the pre- 
vious Section, we proceed to numerically evaluate the full 
Hamiltonian Eq. (13), in the interaction picture, taking 



into account the thermal reservoir for both the atom and 
the field, such that 



P 



(sa) _ 



V,p 



(sa) 



+ ^ (2a,,p(-)a,, - - p^^'^V,,) , (20) 



where 



and Ua and ris are the average number of thermal photons 
for the bosonic field and atomic reservoirs, respectively. 
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FIG. 1: Evolution of tlie numerically computed fidelity for 
both cases of absolute zero (solid line) and finite tempera- 
ture, with the average number of thermal photons n — 0.01 
(dashed line) and n = 0.1 (dotted line). The constant value 
J- = 0.5 indicated by the dashed-dotted line refers to the 
case where the coupling g between the atomic system and the 
cavity mode is null. 

We emphasize that, differently from the master equa- 



tion ( 16 ) which has been derived under the RWA ap- 



proximation leading to the effective Hamiltonian (151, 



the above master equation ( 20 ) describes exactly the dy- 



namics of the whole system thus allowing us to quan- 
tify the errors introduced by our approximations - RWA 
in Eq. (151 and adiabatic elimination of the field vari- 



ables in Eq. (17). Since we are assuming the atomic 
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frequency near the resonance with the cavity field fre- 
quency, both their reservoirs wiU have the same average 
thermal photons, such that from now on we take fia = 
ris- By projecting the above equation in the Fock {\n)} 
and electronic {|e) , \g)} bases, we are lead to an infin- 
ity set of coupled equations for the matrix elements. To 
solve numerically this system of infinity coupled differ- 
ential equations we must truncate the Fock basis some- 
where. The strong decay rate k allows us to safely do it 
since the matrix elements corresponding to highly excited 
Fock states are virtually zero. We then numerically solve 
the master equation ( 20 ) following the method presented 
in [15] . Since we are interested in the evolution of the 
atomic two-level system only, we trace numerically over 
the bosonic field variables. Therefore, after solving the 
full master equation we end up with a density matrix for 
the two-level system with elements 



N 



(21) 



n=0 



where p'^^ = {i\ pat 

Here we analyze the robustness of the protected atomic 
superposition given by Eq. (18). Our strategy con- 



sists in computing the robustness of the protected state 
under thermal effects considering the following different 
regimes: i) <I>(t) = constant, which corresponds to the 
static case, ii) by allowing to slowly vary with time, 
corresponding to Ac/Teng <C 1, and iii) allowing Ac{t) to 
rapidly vary in time, corresponding to Ac/Teng > 1- The 
robustness of the protected state is computed through 
the fidelity T{t) = TT[\+{t)) {+{t)\pat{t)] where patit) 
follows from the numerical solution of the full master 



equation (20). In units of 7, we assumed the reasonable 
decay rates k = 10^7, g = 10^7 and = 2 x 10'^7. 

Starting with the static case, we consider Ac{t) — 
with constant (f>c, such that |-|-) = (|e) -I- e^'^" \g)). As- 
suming (l)c — and the TL atom prepared in the ground 
state, in Fig. 1 we present the numerical results of the 
fidelity against the parameter •yt for both cases of ab- 
solute zero (solid line) and finite temperature, with the 
average number of thermal photons n = 0.01 (dashed 
line) and n = 0.1 (dotted line). As expected, we ver- 
ify that the temperature effects reduce substantially the 
fidelity of the atomic superposition. To stress the effec- 
tiveness of our protocol, we call the attention to the case 
where the coupling g between the atomic system and the 
cavity mode is null. In this case we observe that the fi- 
delity remains J- = 0.5 (dashed-dotted line) showing that 
the cavity mode is a crucial ingredient, together with the 
classical pumping, to protect the desired superposition 
state. 

In Fig. 2 we plot the fidelity of the protected state 
|-|- (t)) undergoing slow and fast evolutions. We have 
considered three distinct values for the detuning between 
the atom and the classical field: Ac = IO7, leading 
to a slow evolution with Ac/T^ng — 0.1 (solid line), 
Ac = 10^7, departing from the adiabatic regime with 
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FIG. 2: Fidelity of the protected state |-|- {t)) undergoing adi- 
abatic and nonadiabatic evolutions at zero temperature. We 
consider the values Ac/Teng — 0.1 (solid line), Ac/Teng — 1 
(dashed line), and finally, Ac/Teng — 3 (dotted line). 



Ac/Teng — 1 (dashed line), and finally Ac = 3 x 10^7, 
with Ac/Tcng — 3 (dotted line). All these curves in Fig. 
2 were plotted considering n = 0. As we are concerned 
with slow and fast evolutions dictated by the time varying 
parameter $(t) = (pc — Act, we focused our attention in 
Ac{t) taking (j)c — 0. We stress that when the condition 
Ac/Tcng ^ 1 is weakened, meaning that some parame- 
ters are rapidly varying in time, we found that although 
the equilibrium is reached more slowly, the fidelity does 
not drop off quickly, attaining a value about 0.9 even for 
Ac/Teng — 3, Corroborating again the effectiveness of our 
protocol. 
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FIG. 3; Fidelity of the protected state |-|- (t)) undergoing 



nonadiabatic evolutions, with Ac = 3 x 10 7 {Ac/Te 



3), 



for absolute zero (solid line) and finite temperature, with the 
average number of thermal photons n = 0.01 (dashed line) 
and n = 0.1 (dotted line). 

In order to see the temperature effects in the nona- 
diabatic evolutions, in Fig. 3 we draw the curves of 
the fidelity for Ac = 3 x 10^7 {Ac/Teng ^ 3) for dif- 



6 



ferent values of mean number of thermal photons, i.e., 
for absolute zero (solid line), n = 0.01 (dashed line) and 
?i = 0.1 (dotted hue). Again, as in the static case (Fig. 
1), we observe that the fidelity decreases with the in- 
crease of the mean number of thermal photons, showing 
that the temperature is the most important source of de- 
coherence in our present protocol. This conclusion is also 
drawn from Fig. 4, where a functional dependence of the 
fidehty of the protected state |+(i)) undergoing nona- 
diabatic evolutions against the reservoirs mean photon 
number is displayed. From this figure, where we have 
used g — IO7 (sohd line), 207(dashed line), 5O7 (dotted) 
and 10^7 (dashed-dotted line), we note that the fidelity 
decays faster with the increasing number of thermal pho- 
tons, as we have mentioned above. 

1i ~ ~ ~ 1 
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20 y 
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0.6- 



0.5' ' ' ' ' 

0.5 j 1.5 2 

n 

FIG. 4; Fidelity of the protected state |-l-(t)) undergoing 
nonadiabatic evolutions, versus the average thermal photons 
n characterizing the finite temperatures with g = IO7 (solid 
line), 207(dashed line), 5O7 (dotted) and 10^7 (dashed-dotted 
line). 



V. CONCLUDING REMARKS 

In this paper we deepened the analysis of the engineer- 
ing reservoir technique we have proposed in Ref. [llj . 
This proposal has been accomplished by deriving a time- 
dependent master equation which leads to a decoherence- 
free evolving superposition state which can be nonadia- 
batically controlled by the system-reservoir parameters. 
In the present contribution we have provided a gen- 
eral recipe to engineer arbitrary effective reservoirs for 
a fermionic system by manipulating its interaction with 
a bosonic mode. 

More specifically, we showed how to protect a superpo- 
sition state of a two-level ion trapped into a leaky cavity. 
The robustness of our scheme was analyzed by means 
of the fidelity considering the case where the protected 
state does not depend on time, as well as the case of 
slowly and rapidly time varying evolutions. To support 
the approximations used to derive our analytical results, 
we have numerically solved the full master equation, ob- 
taining an excellent agreement with the approximations 
we have carried out. In the present contribution we also 
analyzed the temperature effects of the reservoirs on the 
fidelity of the static and time varying protected states. 
We concluded that the temperature is the most impor- 
tant source of decoherence in our present protocol. 

We hope that this contribution can be useful for infor- 
mation processing with trapped ions inside optical cav- 
ities, for example for implementing Deutsch algorithm 
|17( I18j or universal dissipative quantum computing |19j . 

We wish to express our thanks to FAPESP, FAPEMIG, 
CAPES, CNPq, and the Brazihan National Institute 
of Science and Technology for Quantum Information 
(INCT-IQ), for the financial support. 
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